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Quasi-resonant van der Waals interaction between non-identical atoms 
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We present a time-dependent quantum calculation of the van der Waals interaction between a 
pair of dissimilar atoms, one of which is initially excited while the other one is in its ground state. 
For small detuning, the interaction is predominantly mediated at all distances by the exchange of 
doubly resonant photons between the two atoms. We find that it presents both temporal and spatial 
oscillations. Spatially oscillating terms depend on the resonant frequencies of both atoms, while the 
frequency of the time oscillations is given by their detuning. We analyse the physical content of our 
findings and discuss to what extent previous conflicting stationary approaches provide compatible 
results. A proper account of causality is found essential in order to obtain the correct result. 


Dispersion forces between neutral atoms are often in¬ 
terpreted as a result of the quantum fluctuation of both 
the electromagnetic (EM) field and the atomic charges 
0, i- A prominent example of those are van der 
Waals (vdW) forces acting between neutral atoms and 
molecules, which are important in atomic and molec¬ 
ular interferometry where they influence the measured 
interference pattern Q. In Quantum Information, vdW 
forces between Rydberg atoms produce a Rydberg block¬ 
ade which may be exploited to realize quantum gates 
0]. In biophysical and chemical processes vdW forces 
are known to play a crucial role for the stability and as¬ 
sembling of molecules [H. 

At zero temperature, two atoms in their ground states 
undergo a series of virtual transitions to upper levels. It 
is the coupling of the charges of each atom to the quan¬ 
tum EM field that induces the correlation between their 
transient dipole moments, giving rise to a non-vanishing 
vdW interaction. The lifetime of the virtual atomic tran¬ 
sitions is very short in comparison to ordinary observa¬ 
tion times and thus, the use of stationary quantum per¬ 
turbation theory is well justified for the calculation of 
this interaction Q. For short interatomic distances r in 
comparison to the relevant transition wavelengths (non- 
retarded regime) the interaction scales as r~^, while for 
large distances (retarded regime) it goes like r 0il. 

The situation is different for excited atoms. First, ex¬ 
cited states are unstable and present finite lifetimes. This 
implies that, generically, the interaction between excited 
atoms must depend on time. Second, if any of the tran¬ 
sitions from the excited to lower atomic levels is relevant 
to the interaction, the exchange of resonant photons be¬ 
tween the atoms must be considered. The energy of the 
interaction mediated by resonant photons is usually re¬ 
ferred to as resonant van der Waals potential in the lit- 
i 3-13|. In the retarded regime the resonant po¬ 


erature 


tential overtakes by far the non-resonant one. It is in this 
regime that different approaches yield conflicting results 
concerning the spatial oscillations of the interaction [H¬ 


2l|. This long-standing problem is the main motivation 


of the present Letter. 

In the following, we address the time-dependent quan¬ 


tum computation of the interaction between two dissim¬ 
ilar atomic dipoles, one of which is excited. The excited 
atom is taken of type A while the atom in its ground 
state is considered of a different type B. Without loss of 
generality we approximate the atoms by two-level sys¬ 
tems of resonant frequencies uja and ujb respectively, 
with respective linewidths T^i and T^. Further, in or¬ 
der to ensure the perturbative nature of the calculation 
and to avoid resonant energy transfer we set the detuning 
Aab = iOA — OJB such that |Aab| > (F^ +Tb)/‘^ and 
|A^b| 3> l}V{T))/h, with WlT) being the interaction 
Hamiltonian at the time of observation, T. Since the 
observation is made for atom A excited, we must have 
T < 2ttT^ . Lastly, we assume without much loss of 
generality Ta,b < lAyisj <C uja,b, which is easily met by 
pairs of alkali atoms. We will refer to this condition as 
quasi-resonant. We will see that it allows for a great re¬ 
duction in the number of calculations and makes the res¬ 
onant potential dominant at all distances. We will show 
that the interaction energy oscillates both in time and in 
space. It contains time-independent terms which oscil¬ 
late in space with frequency 2u;a/c, and time-dependent 
terms which oscillate in time with frequency Aab and in 
space with frequency 2ujbIc. We compare our results to 
previous conflicting approaches and discuss in detail to 
which extent they provide compatible results. 

We aim at computing the EM energy of atom A due to 
the presence of atom B. To this end we apply standard 
time-dependent quantum perturbative techniques in the 
electric dipole approximation (^ . At any given time T 
the state of the two-atom-vacuum system can be written 
as l'I'(T)) = 1[J(T)|4'(0)), where the state of the system at 
time 0 is |4'(0)) = |A+) <8) |i?-) |0.^). In this expression 

(A, i?)+__ label the upper/lower internal states of the 
atoms A and B respectively and |0.y) is the EM vacuum 
state. U(T) denotes the time evolution operator in the 
Schrbdinger representation, 

U(T) = rexp|-ift-^^ dt 

In this equation Ha Hb is the free Hamiltonian of 
the internal atomic states, fUjJA\A+){A+\-\-hujB\B+){B+\, 


Ha+Hb+ Hem + H }• 
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while the Hamiltonian of the free EM field is Hem = 
/ia;(ak,eajj. + 1/2), where w = ck is the photon 

frequency, and the operators aj^ ^ and ak,e are the cre¬ 
ation and annihilation operators of photons with mo¬ 
mentum Kk. and polarization e respectively. Finally, the 
interaction Hamiltonian reads W = Wa + Wb, with 
Wa. B — B ■ E(Ra,s). In this expression Aa,b are 
the electric dipole operators of each atom and E(Ra,s) 
is the electric field operator evaluated at the position of 
each atom, which can be written in the usual manner as 
a sum over normal modes as @ 


E(Ra,b) = ^E^”)(Ra.b)+e'^+)(Ra,s) 

hck , 


= *E 


2Veo 


eukc 


ik-R/i.B * t„-ikRA 

C U-1,C 


where V is a generic volume and e|/^^ denote the an¬ 
nihilation/creation electric field operators of photons of 
momentum Kk, respectively. While the internal atomic 
and EM degrees of freedom are quantum variables, the 
position vectors Ra.b are classical variables. We empha¬ 
size here that we do not make further simplifications to 
these potentials. In particular, we do not replace the EM 
response of any of the atoms by its ordinary polarizabil¬ 
ity, as it is the case in R ef. [2fl| . 

Next, considering W as a perturbation to the free 
Hamiltonians, the unperturbed time-evolution opera¬ 
tor for atom and free photon states is Uo{t) = 
exp [—ih~^{HA + Hb + HEM)t]- In order to make con¬ 
tact with a realistic setup, we imagine that atom A starts 
being excited at time — r by a laser pulse of duration 
T < |A)^g|. This fixes our temporal resolution and im¬ 
plies that at time ~ 0 the initial state |'I'(0)) is well- 
defined within a time interval of the order of r. We 
are now ready to compute the EM energy of atom A 
due to the presence of atom B at any time T such that 

{Wa{T)) = -(^'(0)|Ut(T)dA • E(RA)U(r)|^'(0)). (1) 


The above expression admits an expansion in pow¬ 
ers of W which can be developed out of the time- 
ordered exponential equation for U(T), U(r) = 
Uo(T)Texp l[jJ(t)WUo(t)dt. At leading order, Eq.([T]) 
contains a series of terms of fourth order in W where 
an electric field operator creates/annihilates a photon 
at time T at the position of atom A. They correspond 
to the twelve well-known time-ordered diagrams of FiglT] 
[1) [l3 • the time-dependent approach, each diagram 
contributes to (Ha(T)) with two terms in which Wa is 
flanked by two U-matrices which make the system evolve, 
in opposite time directions, from the initial state to two 
different states at time T, which differ from one another 
in the state of atom A and in the number of photons by 
one unit. In quasi-resonant conditions, the greatest con- 
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FIG. 1: Diagramatic representation of the twelve time- 
ordered processes which contribute to {Wa{T)) at the lowest 
order in W. The time variable runs along the vertical. 


tribution to {Wa{T)) comes from diagram (a), in which 
two doubly resonant photons are exchanged one after the 
other. Doubly resonant photons are those emitted by one 
of the atoms in its upper level and absorbed by the other 
atom in its lower level, while for non-resonant photons 
the emission/absorpsion processes are inverted. Lastly, 
simply resonant photons are those emitted and absorbed 
by both atoms in either their upper or lower levels. In 
addition, the diagrams (&) — (/) of Fig[Tl which contain 
both doubly resonant and non-resonant photons, provide 
terms which make it possible to extend the frequency in¬ 
tegrals of diagram (a) into the negative domain. Their 
contribution is indeed essential for establishing causal¬ 
ity (cf. Ref. 0 ). All the other contributions from these 
and from the rest of diagrams are at the most of the 
order of Aab/wa.b times smaller and hence negligible. 
Putting everything together, transforming the sums over 
photon momenta into continuum integrals and imposing 
the causality condition T > 2i?/c with R = R^ — Ra, 
we find at leading order, 


(Wa(T)} 


Vk^dk 


Vk'^dk' 
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dD' 


/O 
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i{^{0)\Uo{-T)\'^{0))e{T-2R/c) dt dt' 

Jo Jo 

[ dt"(T'(0)|dA-E^7^(RA)Uo(r-t) (2) 

Jo 


X dB ■ Ej^t)(RB)Uo(t - t')dB ■ E)^-^(Rb) 

X Uo(t' - t")dA ■ E^+)(RA)Uo(t")|^(0))l + [fc o k']l 


A-)t 


The time integrals of the time-evolution operators in 
Eq-© yield a series of terms with poles along the real 
axis. 


_c_ c cos (AabT) 

^ABik - kA){k'- kA) AAB{k - kB){k'- kB) 
cos [(w — UJa)T] cos [{uj' — UJa)T] 

~ {k - kA){k - kB){k'- k) {k' - kA){k' - kB){k' - k)' 

Further, the development of this expression contains 
terms in which both photons resonate either with the 
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transition of atom A or with the transition of atom B 
only. This is a direct consequence of energy conserva¬ 
tion. Upon integration in frequencies, the former terms 
are time-independent while the latter oscillate in time 
as ^ cos AabT. Important is the fact that only the first 
term in Eq. ^ arises in the stationary approach Sim- 
However, the integration in frequencies of the third and 
fourth terms provides additional time-independent con¬ 
tributions which are missing in the stationary approach. 
Lastly, replacing the time integrals in Eq.([2]) with the re¬ 
sult (| 31 ) and integrating in orientations and frequencies, 
we obtain -see Appendix, 


{Wa{T)) = I3P<1 - 

07 y.. 

-b cos{2kAR) + -^P-kA[l3’-^ 

- k\R‘^a^^pP‘^] sin (4) 

- - klR^{P^^pP^ + 2a^^pP’^) 

o/y.. 

+ k%R‘^a^^aP^] cos (2fcsi? + AabT) - 
X [/3*J/3P9 - klR^a^^pP’^] sin {2kBR + AabT) 

H-[1 + + {kA,BR)^\0 {AAB/^a,b) + 


where = /if/i^/r^/r^/[(47reo)^/iAAB], = 

,B 


(A-|d. 4 |A+), /i^ = {B^\dB\B+) and = <5*^' - 
= 5^^ — R^R^/R^. It is worth stressing 
that the Heaviside function in Eq.([2|) together with the 
time-order prescription do not only guarantee causality, 
but also determine univocally the contours of integra¬ 
tion over frequencies in the complex plane when taking 
the principal value (see Appendix). The last term in 
Eq.® indicates the order of the leading corrections to 
the dominant doubly resonant photon exchange terms of 
Eq.(l2|) 22|. As anticipated, the time-independent terms 
of Eq.Q oscillate only in space with frequency 2kA- On 
the contrary, the time-dependent terms oscillate in time 
with frequency Aab and in space with frequency 2kB- 
Only for large integration times, 5T » their time 

average vanishes. In the short time limit, T -A 2R/c, 
(Wa{T)) vanishes identically at our order of approxima¬ 
tion. This is a consequence of the fact that, in order to 
establish the interaction, it is necessary that the excita¬ 
tion be transferred actually to atom B. For T > R/c, 
the probability of excitation of atom B oscillates in time 
as K«'(T)|A_)®|H+)®|0^)|2 ^sin^ [AAB(i?/c - T)]/2], 
being maximum for the first time at T = RjcA 7r|A)(g|. 
Correspondingly, (WAiT)) becomes maximum for the 
first time at T = 2i?/c-|- 7r|A)^g|. The lapse R/c be¬ 
tween these two times is the time for a photon to travel 
back from to R^i after the excitation of atom B. 

A long-standing debate exists in the literature concern¬ 
ing the spatial oscillations of the two-atom interaction in 
the retarded regime when one of the atoms is excited 


3“21|. The existence of spatial oscillations is indeed 


supported by experiments [2f 


According to our 


findings, for kA.BR ^ 1 and T > 2i?/c, the interaction 
oscillates both in time and in space as 


{Wa{T)) ~ cos {2k aR) (5) 

- k% cos {2kBR + AabT)] 

- —^a^^aP'^k\ sin [Aab{R/c - T/2)] 

X sin [kA{R + cT/2) + kB{R — cT/2)]. 


From the last expression we read that, at fixed time, the 
interaction is modulated by long-range oscillations of fre¬ 
quency Aab /c, while short-range oscillations take place 
at frequency kA + kB- Also as a function of time the in¬ 
teraction is modulated by oscillations of frequency Aab ■ 
In Fig 12] we plot the energy of the interaction between 
two alkali atoms, one of ®^Rb which is excited to the 
state 5 Pi/ 2 and another one of which is in its ground 
state, in the retarded regime. 

In contrast to our result, the stationary approach of 
Power and Thirunamachandran in Ref.[l4| predicts no 
oscillations for {Wa) in the far field. The key point in 
their calculation is the addition of small imaginary parts 
to the resonant frequency of atom A in such a way that 
poles get shifted off the real axis. They used the pre¬ 
scription that a positive/negative imaginary part must be 
added for emitted/ absorbed photons in order to account 
for the finite linewidth of the excited atom. In particular, 
for Aab 'C ia)a,b the dominant term in their stationary 
calculation is the first one in Eq.(l3|), but with the real 
poles shifted as [AAB{k — kA— iT]/c){k' — kA + iri/c)]~^, 
7/ —>■ 0^. After integrating in orientations an analogous 
equation to Eq.Q [6|, they must have obtained for the 
energy in the far field limit, ;j> 1, 

7/ r+oo r+oo 

-/ dk dk'Pk'^ (6) 

47r2R2 J_^ J_^ 

^i{k-\-k^)R _j_ ^—i{k-\-k')R ^i{k—k')R ^—i{k — k^)R 

_ 

{k — kA — irj/c){k' — kA + itl/c) ’ 


with ry ^ O’*". Since the pole in k lies on the upper half of 
the complex plane and the pole in fc' lies on the lower half, 
the only nonvanishing contribution to the above integral 
comes from the term proportional to As the 

real part of the poles is in both cases fc^, taking the 
limit ?7 —> O'*" the exponent vanishes after evaluating the 
residues and the integral yields the non-oscillating result 

In the previous stationary calculation of McLone and 
Power [13 and Gomberoff et al. the poles in Eq.® 
were not shifted. As a result, when taking the principal 
value of the integrals in Eq.® with 77 = 0, the four expo¬ 
nentials in the numerator contribute as ^ 2 cos 2k aR + 
2 = 4cos^ kAR after adding up the residues, yielding the 
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FIG. 2: Graphical representation of the interaction between 
a ®^Rb atom in state 5Pi/2 {kA = 2n 12578.95 cm“^) and 
a atom in its ground state, for kAR 1. The black 
line corresponds to a snapshot of the interaction at time 
3.0 ■ s, where the contributions of the D1 and D2 tran¬ 

sition lines of the atom {kg — 27r 12985.17 cm~^ and 
k% = 27r 13045.876 cm“^ respectively) add up approximately 
in-phase. A long-range period ctt/IAasI, with AabIc = 
kA — (kg -I- k%)/2, is identified. The red line corresponds to 
the time-independent result of the causal-adiabatic approxi¬ 
mation. Average over dipole orientations has been taken. 


oscillating result cos^ k^R- As mentioned 

after Eq. some time-independent terms are missing 
in the stationary calculation, which explains the discrep¬ 
ancy of this result with the time-independent component 
of ours in Eq.®. In fact, the correct time-independent 
result can be obtained by switching on adibatically the in¬ 
teraction potential, W, within the time-dependent causal 
approach. That is, by replacing WA,B(t) in Eq.® with 
?7 —S' 0 +, and extending the lower limit of the 
time integrals to —oo-see Appendix. This results in an 
equation analogous to Eq.® but with both poles shifted 
along the imaginary axis in the same direction an amount 
iTq/c, which yields k\a'‘^aP'i cos 2kAR for » 1 . 

Recently, Safari and Karimpour have published a let¬ 
ter [l^ where they claim to obtain for UaR ^ 1 the same 
oscillating behaviour as Gomberoff et al. 0. However, 
a straightforward comparison of Eq.(19) of Ref. 0 and 
Eqs.(14,26) of Ref. [1^ reveals that this is indeed not 
the case. Whereas the result of the latter is the one 
outlined above, cos^ kAR, the authors of the former 
have found ^ cos2k aR, despite the fact that both ap¬ 
proaches are based on fourth order stationary perturba¬ 
tion theory. The origin of the discrepancy is in the alge- 
braic manipulation inherited by the authors of Ref.[l9| 
from Ref.Q. In the Appendix B of Ref.j^ the authors 
have tried to express the total contribution of the twelve 
diagrams of FigU] as a single frequency integral whose 
integrand is a function of the ordinary polarizabilities 
of the two atoms. In doing so by means of Eq.(B2) 
of Ref.0, the authors have replaced effectively the de¬ 
nominator of Eq.®, which is a symmetric and separa¬ 
ble function of k and k' for 77 = 0 , by the expression 
[Aab(A: — kA)]~^[l/{k' — k) + l/{k' + fc)], which is nei¬ 


ther symmetric nor separable. As a consequence, that 
replacement makes the frequency integrals depend arbi¬ 
trarily on the order of integration. Next, integrating in k' 
first and in k later, one obtains k\ cos 2k aR, 


which agrees with Eq.(19) of Ref.|19j| for kAR 3> 1, 
Aab ^ ^A,B, upon averaging in atomic orientations. In¬ 
terestingly, this result equals the time-independent term 
of Eq.®. However, this coincidence can only be ac¬ 
cidental, since the above replacement and the subse¬ 
quent prescription on the order of integration are neither 
connected to the time-dependent terms of Eq.® which 
cause the actual discrepancy with respect to the result of 
Refs. 15, 1^ nor to the causal-adiabatic approximation. 

It is worth noting that while we have invoked the exis¬ 
tence of finite lifetimes ~ ^^b order to impose phys¬ 
ical constraints on the detuning Aab and on the obser¬ 
vation time T, no explicit reference to these quantities 
appear in our expression for (ITa(T)). As a matter of 
fact, only the emission through the exchange of reso¬ 
nant photons between the two atoms has been implic¬ 
itly accounted for in our calculation of {Wa{T)). How¬ 
ever, our calculation lacks the inclusion of the spon¬ 
taneous emission of each atom into free space, whose 
rates are k^ g/^neoh, respectively. The processes 
corresponding to the latter phenomenon are generally 
unimportant in comparison to those depicted in FiglT] 
since their leading contribution to {Wa{T)) is of order 
0{W^) ^ MaMb; f^Ak-B FigO They might only be 
relevant for the case that the lifetimes are of the order of 
the temporal frequency of the interaction, Ta.b | Aab|, 
but they cannot affect in any case the oscillatory be¬ 
haviour found here for the terms of order This 

argument opposes to the reasons given in Ref. [1^ to add 
imaginary shifts to the real poles at 0{W'^). We final¬ 
ize by mentioning that Berman has shown in a recent 
publication how to introuduce spontaneous emis¬ 
sion in a consistent manner in an adiabatic approxima¬ 
tion. Assuming that atom A is excited adiabatically 
with t,5T ^ he has obtained the correct time- 

independent result. 

In this Letter we have shown that the van der Waals 
interaction between two dissimilar atoms, one of which is 
initially prepared in an excited state, presents generically 
oscillations both in time and in space. In quasi-resonant 
conditions the interaction is dominated at all distances 
by the exchange of doubly resonant photons between the 
two atoms. It is modulated in space by long-range oscilla¬ 
tions of frequency Aab /c, while short-range oscillations 
take place at frequency kA + ks- The time frequency 
is Aab, which determines the rate at which the excita¬ 
tion is transferred to atom B. In the retarded regime the 
interaction takes the form of Eq. ® . Only for large inte¬ 
gration times, ST :$> |A(^g|, that expression reduces to a 
time-independent term which oscillates in space with fre¬ 
quency 2kA- The latter is also equivalent to the result of 
the causal-adiabatic approximations (see Appendix and 
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Ref.[2l|). It does not agree, however, with the result of 
stationary perturbation theory 15l. Il6|. 



FIG. 3: Dominant diagrams at order 0{W^) which incorpo¬ 
rate in the calculation of {Wa{T)) the effect of photon emis¬ 
sion into free space. In (i) and {ii) the photon 7 is emitted 
into free space from atom A, whereas in {iii) it is emitted 
from atom B. 

We thank Paul Berman and Marie-Pascale Gorza for 
fruitful discussions. Financial support from ANR-IO- 
IDEX-0001-02-PSL and ANR-13-BS04-0003-02 is grate¬ 
fully acknowledged. 


Appendix: Integration of Eq. ([^ in the complex 
plane and the causal-adiabatic approximation 


In the first place, we note that all the time-evolution 
free-propagators in Eq.(2) are retarded, Uo(t“ — t^) with 
= 0,t",t',t,T. This is so because the time-order 
prescription on the time integrals implies T > t' > t" > 
0 , hence holds in all the Uo(t“ — t^) oparators of 

Eq.(2). It is the time-ordering together with the Heavi¬ 
side function 0(T —2i?/c) in the integrand of Eq.@ that 
guarantee causality. 

Substituting in Eq.Q the expressions of the operators 
E(±)(Ra, b) as sums over normal modes, and making use 
of the identities @ X]e(k) “ kikj/k'^ and 


/*47r 


an, 

47r 


{Sij — kikj/k‘‘)e^^'^ = aij sinkR/kR 


■ fiij [cos kR/(kR)^ — sin kR/(kR)" 


(7) 


that we will denote by J^ij{kR), we can write Eq.Q as 


I*+oo 


(WAiT)) / dkk^r^ikR) 

oTT CqC J_qq 

p-\-CXD nT 

X / dk'k'^RP’^ik'R) / dt0(T-2R/c) 

J — 00 Jo 


' —00 


dt' 


dt" 


(* 


^iu]AT -i{T-t)u 


( 8 ) 




where the integrand of the time-ordered integrals is the 
product of the retarded time-evolution free-propagators, 
a = e^/ 47 reofic is the fine-structure constant and e is the 


electron charge. Next, integrating in time we end up with 

/ + 00 

dkk^r^ikR) 

-00 

/ +00 

dk' k'^RP‘^{k'R)e{T - 2R/c) 

-00 

r c /^AB clA ab COs(Aab^) 

L(k - kA)(k'- kA) (k - kB)(k'- ks) 
cos [(w — a;^)r] 


(k — kA)(k — kB){k' — k) 
cos [{uj' — uja)T] 

{k'- kA){k' - kB){k' - k)y 


where straightforward application of perturbation theory 
enforces to taking the principal value (P.V.) of the inte¬ 
gral in Eq.®, as intermediate states with energies equal 
to that of the initial state must be removed from the 
sums (integrals in the continuum) in order to avoid zeros 
in the denominators. 




fko—’n/c 




k — ko 


Hk) 

k — ko 




dk 


Hk) 

fko-\-r}/c ^ ^0 


dk 


^+00 


= P.V. 


R{k) 


k — h 


0 


dk as ry —> O'*'. (10) 


Having impossed causality by time-ordering the integrals 
of Eq.® and multiplying the integrand by 0(T — 2i?/c), 
there is no need to add arbitrarily small imaginary parts 
±i0+ to the real poles of the integrand in Eq.(l9|) in or¬ 
der to ensure that the propagators are retarded. Also, 
we have argued in the Letter that spontaneous emission 
enters at higher order, hence it cannot provide imagi¬ 
nary shifts ±iTA,B/c to the real poles of Eq.® either, 
in contrast to Refs.[l3, [Uj. Lastly, the contours of inte¬ 
gration for each of the terms of the frequency integrals in 
the complex plane is univocally determined by the con¬ 
dition T > 2Rlc. Each term contains complex exponen¬ 
tial factors of k and k' of the form h{R,T)^ 

with fi^ 2 {R,T) being generic linear functions of R and 
T. That implies that the integration contours must be 
closed by infinitely large semi-circles either in the upper 
half plane for fi^ 2 {R,T) > 0 , or in the lower half plane 
for /i, 2 ( 7 ?, T) < 0, which is ultimately determined by the 
condition T > 2i?/c. It is worth mentioning that the in¬ 
tegrand of Eq. ( 0 ) is invariant under the exchange fc -O- fc', 
which implies that the total integral is independent of the 
order of integration, as expected [i^. Thus, (IE 4 (T)) is 
unambiguously given by Eq. ([9]) . The final result is given 
in the Eq.(4) of the Letter. 

Finally, the causal-adiabatic approximation referred to 
in the conclusions of the Letter corresponds to adding 
up a factor +* '> to the integrand of Eq.®, with 

77 —>■ 0 +, and to extending the lower limits of the time 
integrals to — 00 . After performing the time integrals one 
obtains the T-independent terms of Eq.® (excluding the 
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Heaviside function) with a small imaginary part —irj/c 
added to both real poles, 


OL ABBA 

47r3eoe2^* , 

X JP^{k'R) 


(' + 00 


dfc k^r^{kR) 

C 


r*+oo 


^Asik -kA- ir]/c){k' - kA - ir]/c) 


dk'k'^ 

■ ( 11 ) 


We also note that Eq.(6) in the Letter corresponds to the 
T-independent terms of Eq.® of which only the terms 
proportional to and aP‘^ have been retained in the 
J^-functions, and in which small imaginary parts ±ir]/c 
have been added to each real pole respectively. 
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